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Dear  Tatjana  and  Erin, 

I’m  writing  because  I  am  a  bit  confused  about  when  my  annual  report  for  Award  No.  FA9550-08- 
1-0069  is  due.  A  document  I  was  sent  when  I  received  the  award  seems  to  indicate  that 
nothing  is  due  until  Feb.  2009  (see  the  end  of  attached  document  123082_2.tif).  However,  my 
recollection  from  previous  years  is  that  annual  reports  for  AFOSR  have  been  due  in  the  Fall. 
Please  let  me  know  if  you  need  something  before  Feb.  28  2009.  We  have  been  getting  great 
results  in  the  lab  and  it  would  be  easy  to  put  something  together  immediately  if  you  need  it. 

I  also  wanted  to  take  this  opportunity  to  ask  whether  everything  looks  good  for  funding 
Option  1  of  this  award.  Please  let  me  know  if  there  are  any  anticipated  problems  or  delays. 

I  also  wanted  to  give  a  brief  update  on  the  results  we've  been  getting  in  the  laboratory. . . 

We  have  been  experimentally  investigating  three-state  mixtures  of  degenerate  Fermi  gases. 
Using  these  mixtures  we  can  confirm  theoretical  predictions  for  the  quantum  3-body  problem. 

In  the  longer  term  (~1  to  2 

yrs.)  we  will  confine  this  mixture  in  an  optical  lattice  and  simulate  the  SU(3)  Hubbard  model 
and  high  density  quark  matter  described  by  quantum  chromodynamics. 

We  have  been  investigating  the  stability  of  the  3-state  Fermi  gas  with  respect  to  decay  due 
to  three-body  recombination.  The  three-body  recombination  rate  is  an  observable  we  can  use 
to  compare  to  theoretical  predications  based  on  solutions  to  the  quantum  3-body  problem. 
Understanding  the  stability  of  the  gas  also  has  important  implications  for  the  feasibility  of 
future  studies  of  quantum  chromodynamics  in  this  system.  We  have  been  working  with  Eric 
Braaten,  a  high  energy  physics  theorist  at  Ohio  State  University,  to  interpret  our  data. 

Eric  has  been  able  to  fit  our  results  using  a  three-body  model  with  only  2  free  parameters. 
His  results  indicate  that  our  results  are  in  agreement  with  predictions  made  by  Vitaly 
Efimov,  a  Russian  nuclear  physicist,  over  35  years  ago.  Futhermore,  Eric's  results  predict 
that  we  should  be  able  to  perform  even  more  stringent  tests  of  Efimov's  predictions  in 
experiments  we  are  performing  right  now.  These  results  will  be  of  interest  in  several 
subfields  of  physics  including  nuclear,  high-energy,  atomic,  molecular  and  particle  physics 
as  well  as  quantum  chemistry. 

Our  results  also  indicate  that  the  observation  of  BCS  superfluidity  in  this  system  should  be 
possible.  The  superfluid  phase  in  this  system  would  be  analogous  to  the  color 
superconducting  phase  predicted  to  occur  inside  neutron  stars.  The  possiblity  for  such 
future  studies  was  highlighted  in  the  October  issue  of  Physics  Today  (see  attached  article). 
We  have  submitted  our  recent  results  to  PRL  (see  attachment). 

We  also  recently  submitted  our  proposal  for  cooling  fermions  in  an  optical  lattice  to  PRL 
(see  attachment). 

Our  experimental  capabilities  have  improved.  In  addition  to  being  able  to  produce  degenerate 
Fermi  gases  composed  of  2  or  3  spin  states,  we  can  also  produce  molecular  Bose-Einstein 
condensates  in  our  laboratory  which  are  composed  of  Li  dimer  molecules.  This  allows  us  to 
get  to  extremely  low  temperatures. 


l 


We  have  also  developed  and  demonstrated  a  novel  ring  laser  sour^?.  In  collaboration  with 
Stephen  Gensemer,  a  former  postdoc  in  Kurt  Gibble's  lab,  we  have  demonstrated  a  new  method 
for  obtaining  uni-directional  single-frequency  operation  of  a  high-power  ring  laser  by  a 
technique  we  call  self  injection  locking.  This  technique  can  be  applied  to  any  ring  laser 
and  therefore  will  be  of  broad  interest  to  any  research/industrial  application  which  uses 
high-power  single  frequency  ring  lasers.  We  are  in  the  process  of  writing  these  results  up 
and  plan  to  submit  a  paper  to  Nature  Photonics. 

Please  let  me  know  if  you  need  an  official  performance  report  prior  to  Feb  28,  2009.  As  I 
said,  I  can  get  something  to  you  quickly  if  necessary.  If  the  report  is  already  delinquent, 
I'm  sorry  for  the  confusion  and  I  will  get  it  in  immediately. 

Best  regards, 

Ken 


Kenneth  M.  O'Hara 

Assistant  Professor  of  Physics 

Penn  State  University 

104  Davey  Lab,  PMB  195 

University  Park  PA  16802 

Phone:  814-865-7259 

FAX:  814-865-3604 

e-mail:  kohara|3phys . psu . edu 
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ultracold  gases 
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The  unprecedented  control  over  the  interactions  and  pairing  of  ultracold  fermionic  atoms  provides 
insight  into  exotic  strongly  correlated  phenomena  and  illuminates  the  physics  of  superfluidity  in 
metals,  nuclei,  and  neutron  stars. 


Carlos  Sd  de  Melo  is  a  professor  of  physics  at  the  Georgia  Institute  of  Technology  in  Atlanta. 


If  you  ask  an  undergraduate  student,  “What  do 
neutron  stars,  metals,  nuclei,  and  atoms  have  in  common?" 
you  might  hear  the  answer,  "They  are  made  of  neutrons,  pro¬ 
tons,  quarks,  and  electrons."  While  it  is  common  for  students 
to  have  heard  about  the  existence  of  those  microscopic  parti¬ 
cles  and  their  individual  properties  like  charge,  spin,  and 
color,  it  is  far  less  common  to  hear  from  them  that  neutrons, 
protons,  quarks,  electrons,  and  atoms  may  exhibit  a  collec¬ 
tive  and  macroscopic  property  called  superfluidity,  in  which 
a  large  number  of  particles  can  flow  coherently  without  any 
friction  or  dissipation  of  heat. 

Neutron  stars,  metals,  nuclei,  and  ultracold  atoms  have 
revealed  their  amazing  superfluid  state  in  ingenious  exper¬ 
iments  (see  figure  1  for  estimated  critical  temperatures  of 
fermion  superfluids).  The  first  discovery  of  superfluid  be¬ 
havior  was  made  nearly  100  years  ago,  when  Heike  Kamer- 
lingh  Onnes  in  1911  cooled  a  metallic  sample  of  mercury  to 
temperatures  below  4.2  K,  using  liquid  helium-4  as  the  re¬ 
frigerant,  and  found  to  his  astonishment  that  the  sample 
conducted  electricity  without  dissipation.  That  dramatic 
drop  in  electrical  resistivity  to  an  essentially  zero-resistance 
state  was  coined  superconductivity,  also  known  as  charged 
superfluidity. 

Revealed  through  Kamerlingh  Onnes's  push  to  reach 
very  low  temperatures,  that  incredible  quantum  phenome¬ 
non  emerged  several  years  before  the  formulation  of 
Schrodinger's  equation  and  the  basic  developments  of  quan¬ 
tum  theory.  It  took  nearly  50  years  after  the  discovery  of  su¬ 
perconductivity  for  the  development  of  a  microscopic  theory 
of  the  phenomenon.  In  the  intervening  years,  many  new  su¬ 
perconductors  (charged  superfluids)  were  discovered,  but 
only  one  neutral  superfluid:  liquid  4He. 

Liquid  helium-4:  A  boson  superfluid 

Although  individual  neutrons,  protons,  quarks,  and  elec¬ 
trons  are  all  fermions  with  spin  1/2,  individual  atoms  and 
nuclei  can  be  either  bosons  with  an  integer  total  spin  (0,  1, 
2,  . . .)  or  fermions  with  half-integer  spin  (1/2,  3/2,  5/2,  . . .), 
depending  on  the  number  of  their  constituent  fermions.  The 
neutral  atoms  of  4He,  comprising  an  even  number  of  fermi¬ 
ons  (two  protons,  two  neutrons,  and  two  electrons),  are 
bosons  with  total  spin  0.  An  unusual  maximum  in  the  den¬ 
sity  of  liquid  4He  had  been  observed  as  early  as  1910  (and 


measured  again  in  1922)  in  Kamerlingh  Onnes's  laboratory, 
and  a  possible  discontinuity  in  the  latent  heat  near  the  same 
temperature  was  noted  in  1923  in  the  same  laboratory,  but 
the  significance  of  those  anomalies  was  not  fully  realized 
until  the  late  1930s.  Experimental  evidence  for  superfluid  be¬ 
havior  in  4He  emerged  only  in  1938  when  two  independent 
measurements  of  viscosity  were  reported,  one  by  Pyotr 
Kapitsa  and  the  other  by  John  Allen  and  Donald  Misener; 
the  measurements  indicated  exceedingly  low,  essentially 
vanishing  viscosity  below  2.2  K,  the  temperature  at  which 
the  density  and  latent-heat  anomalies  of  liquid  4He  ap¬ 
peared.  Not  long  after  that  discovery,  Fritz  London* 1  sug¬ 
gested  the  connection  between  superfluidity  and  Bose- 
Einstein  condensation  (BEC). 

The  idea  of  boson  condensation  was  put  forth  by  Albert 
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Figure  1.  Fermionic  superfluidity  is  a  collective  phenomenon 
found  in  neutron  stars,  nuclear  matter,  superconductors,  liquid 
helium-3,  and  ultracold  atoms,  though  they  exhibit  large  dif¬ 
ferences  in  the  critical  temperature  for  the  onset  of  super¬ 
fluidity,  as  shown  here.  The  key  to  their  common  behavior  is 
the  pairing  of  fermions  to  produce  highly  correlated  fluids, 
which  reveal  themselves  in  several  ways:  They  flow  without 
dissipation,  may  exhibit  an  energy  gap  in  the  elementary 
excitation  spectrum,  have  a  nonclassical  moment  of  inertia, 
and  may  exhibit  vortex  lattices  upon  rotation. 
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Figure  2.  Phase  diagram  for  fermionic  superfluids.  To 
see  the  superfluid  Bardeen-Cooper-Schrieffer  or  Bose- 
Einstein  condensation  phases  for  fermions,  it  is  necessary 
to  cool  them  to  a  small  fraction  of  their  Fermi  temperature 
Tf.  When  fermions  pair  to  produce  a  superfluid  state  (yel¬ 
low  region),  they  exhibit  two  characteristic  temperatures: 
the  pairing  temperature  T  „  (pink)  at  which  fermion  pairs 
first  form  and  the  critical  temperature  Tc  (blue)  at  which 
phase  coherence  between  the  pairs  is  established.  In  the 
BCS  regime  of  large  pairs  and  weak  attraction,  the  two 
temperature  scales  are  essentially  the  same.  However,  in 
the  BEC  regime  of  small  pairs  (diatomic  molecules)  and 
strong  attraction,  the  two  temperature  scales  are  very  dif¬ 
ferent.  The  unitarity  limit  occurs  when  the  scattering  pa¬ 
rameter  1  / kFa1  goes  to  zero,  where  kF  is  the  Fermi 
wavenumber  and  as  is  the  scattering  length,  but  the  de 
facto  separation  between  the  BCS  and  BEC  regimes  is 
where  the  chemical  potential  /u  goes  to  zero.  (Adapted 
from  ref.  7.) 
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Einstein  in  1924-25,  after  he  had  read  and  translated  a  paper 
by  Satyendra  Nath  Bose  on  the  statistics  of  photons  and  had 
extended  Bose's  idea  to  the  realm  of  massive  bosons.2  Einstein 
noted  that  a  finite  fraction  (called  the  condensate  fraction)  of 
the  total  number  of  massive  bosons  would  macroscopically 
occupy  the  lowest-energy  (zero  momentum)  single-particle 
state  at  sufficiently  low  temperatures.  For  noninteracting 
bosons,  the  fraction  of  zero-momentum  particles  is  100%  at 
zero  temperature,  but  interactions  in  liquid  4He  are  suffi¬ 
ciently  strong  to  reduce  the  condensate  fraction  to  10%,  even 
at  zero  temperature. 

Some  microscopic  understanding  was  developed  by 
Nikolai  Bogoliubov,  who  demonstrated  that  the  excitation 
spectrum  of  a  weakly  interacting  Bose  gas  is  linear  in  the  mo¬ 
mentum  of  the  excitation  and  that  the  critical  velocity— the 
flow  speed  at  which  superfluidity  is  destroyed  — is  finite.  Bo¬ 
goliubov  concluded  that  weak  repulsive  interactions  do  not 
destroy  the  Bose-Einstein  condensate,  and  that  an  ideal  Bose 
gas  in  its  BEC  phase  has  a  vanishing  critical  velocity  and  thus 
is  not  a  superfluid.  Therefore,  quite  generally,  the  formation 
of  a  Bose-Einstein  condensate  does  not  guarantee  superflu¬ 
idity,  which  is  associated  with  the  existence  of  currents  that 
flow  without  dissipation  and  requires  correlations— that  is, 
interactions— between  bosons. 

Superconductivity  and  the  BCS  theory 

Drawing  from  the  lessons  of  superfluid  4He  and  the  connec¬ 
tion  to  BEC,  Max  Schafroth  proposed  in  1954  that  super¬ 
conductivity  in  metals  was  due  to  the  existence  of  a  charged 
Bose  gas  of  two-electron  bound  states— local  fermion  pairs— 
that  condense  below  a  critical  temperature.  However,  ex¬ 
periments  did  not  seem  to  support  that  simple  picture.  It  was 
not  until  1956  that  a  key  idea  emerged:  Leon  Cooper  dis¬ 
covered  that  an  arbitrarily  small  attractive  interaction  be¬ 
tween  two  fermions  (electrons)  of  opposite  spins  and  oppo¬ 
site  momenta  in  the  presence  of  many  others  could  lead  to 
the  formation  of  bound  pairs.  Since  the  pairing  occurs  in 
momentum  space  and  the  attraction  between  electrons  is 
weak,  the  bound  pairs,  now  known  as  Cooper  pairs,  are 
quite  large  and  thus  different  from  the  local  pairs  envisioned 
by  Schafroth.  The  origin  of  the  glue  between  fermions  was 
argued  to  be  electron-phonon  interactions,  which  could  pro¬ 
duce  an  effective  attractive  interaction  between  electrons 


that  in  turn  could  overcome  their  Coulomb  repulsion. 

The  existence  of  a  single  fermion  pair,  however,  was  not 
sufficient  to  describe  the  macroscopic  behavior  of  supercon¬ 
ductors.  It  was  still  necessary  to  invent  a  collective  and  cor¬ 
related  state  in  which  many  fermion  pairs  acting  together 
could  produce  a  zero-resistance  state.  The  invention  of  such 
a  special  state  of  matter  occurred  in  1957,  when  John  Bardeen, 
Cooper,  and  Robert  Schrieffer  (BCS)  proposed  a  many- 
particle  wavefunetion  corresponding  to  largely  overlapping 
fermion  pairs  with  zero  center-of-mass  momentum,  zero  an¬ 
gular  momentum  (s-wave),  and  zero  total  spin  (singlet).3  As 
emphasized  by  BCS,  their  theory  did  not  describe  the  picture 
proposed  by  Schafroth,  Stuart  Butler,  and  John  Blatt,  also  in 
1957,  in  which  Bose  molecules— local  pairs  of  electrons  with 
opposite  spins— form  an  interacting  charged  Bose  gas  that 
condenses  and  becomes  a  superconductor.4 

One  of  the  most  fundamental  features  of  the  BCS  state  is 
the  existence  of  correlations  between  fermion  pairs,  which 
lead  to  an  order  parameter  for  the  superconducting  state.  In 
the  original  BCS  work,  the  order  parameter  was  found  to  be 
directly  related  to  the  energy  gap  EK  in  the  elementary  exci¬ 
tation  spectrum.  Since  all  relevant  fermions  participating  in 
the  ground  state  of  an  s-wave  superconductor  are  paired,  cre¬ 
ating  a  single  fermion  excitation  requires  breaking  a  Cooper 
pair,  but  that  costs  energy.  Thus  the  contribution  of  elemen¬ 
tary  excitations  to  the  specific  heat  and  other  thermodynamic 
properties  shows  exponential  behavior  ~  exp(-EK/T)  at  low 
temperature  T.  Unlike  the  alternative  theory  of  Bose  mole¬ 
cules,  the  BCS  theory  was  an  immediate  success,  since  it 
could  explain  many  experimental  results  of  the  time  at  a 
quantitative  level. 

After  its  great  success  in  describing  superconductors 
(charged  superfluids),  the  BCS  theory  was  quickly  general¬ 
ized:  to  higher  angular  momentum  pairing;  to  neutral  su¬ 
perfluids  such  as  liquid  3He  (a  fermion  isotope  of  helium  with 
two  protons,  one  neutron,  and  two  electrons),  where  spin- 
triplet  p-wave  superfluidity  was  found  experimentally  in  the 
1970s;  and  to  pairing  mechanisms  beyond  phonon  media¬ 
tion.  But  equally  important  were  applications  of  the  BCS  the¬ 
ory  to  describe  superfluid  phases  of  fermions  in  nuclear  mat¬ 
ter  and  neutron  stars,  where  nuclear  forces  provide  the  glue 
for  fermion  pairs.  Glitches  in  the  rotational  periods  of  neu¬ 
tron  stars  have  been  attributed  to  the  depinning  of  vortices— 
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Feshbach  resonances 


The  evolution  from  a  Bardeen-Cooper-Schrieffer  superfluid  to  a  Bose-Einstein 
condensation  superfluid  cannot  be  studied  in  neutron  stars,  nuclear  matter, 
superconductors,  or  liquid  helium-3,  but  in  ultracold  atoms  it  can  be.  Feshbach 
resonances  are  the  tools  that  allow  the  interactions  between  atoms  to  be  changed 
as  a  function  of  the  applied  magnetic  field.  The  underlying  requirement,  shown 
schematically  on  the  left,  is  that  at  zero  magnetic  field,  the  interatomic  poten¬ 
tials  of  two  atoms  in  their  ground  state  (the  so-called  open  channel)  and  in  an 
excited  state  (the  closed  channel)  be  not  too  different  in  energy.  The  resonance, 
characterized  by  a  divergence  in  the  scattering  length  a,,  occurs  when  the  energy 
difference  AE  between  a  bound  state  with  energy  Em  in  the  closed  channel  and 
the  asymptotic,  threshold  energy  E of  scattering  states  in  the  open  channel 
is  brought  to  zero  by  an  applied  external  magnetic  field  B0.  For  magnetic 
fields  B  close  to  B0,  the  background  scattering  length  is  renormalized  to 
as  «  citjl  +  A B/(B  -  B0)],  where  AB  is  the  width  of  the  Feshbach  resonance. 
When  a^  is  negative,  as  is  positive  close  to  the  resonance  for  magnetic  fields 
smaller  than  B0,  as  seen  on  the  right,  and  two  fermions  can  form  Feshbach  mole¬ 
cules  of  characteristic  size  a,.  The  most  commonly  studied  Feshbach  resonances 
for  s-wave  scattering  occur  at  B0  =  83.4  mT  (834  gauss)  for  lithium-6  and  at 
B0  =  22.4  mT  (224  gauss)  for  potassium-40,  both  stable  fermionic  isotopes. 


fundamental  excitations  of  rotating  superfluids— in  the  stars' 
solid  crust.  For  some  nuclei,  the  energy  gap  in  the  elemen¬ 
tary  excitation  spectrum  and  the  nonclassical  moment  of  in¬ 
ertia  indicate  the  existence  of  a  superfluid  state. 

From  BCS  to  BEC  superfluids 

Although  the  BCS  theory  and  its  generalizations  have  found 
applications  in  several  areas  of  physics  where  the  mecha¬ 
nisms  for  fermion  pairing  are  quite  different,  it  is  intrinsically 
a  weak-attraction  theory.  Nature  was  very  kind  to  the 
BCS  theory  because  it  saved  some  of  its  most  precious 
Fermi  superfluids  to  be  discovered  only  in  the  mid-1980s 
(high-T  cuprate  superconductors)  and  mid-2000s  (ultracold 
lithium-6  and  potassium-40  fermionic  atoms),  when  strong 
deviations  from  BCS  behavior  were  found.  But  prior  to  those 
experimental  discoveries,  a  generalization  of  the  BCS  theory 
was  slowly  developed  to  encompass  the  strong-attraction 
regime  in  which  fermion  pairs  become  tightly  bound  diatomic 
Bose  molecules  and  undergo  Bose-Einstein  condensation. 
The  key  question  addressed  theoretically  was.  Are  the  BCS 
and  BEC  theories  the  endpoints  of  a  more  general  theory  that 
connects  Fermi  superfluids  to  molecular  Bose  superfluids? 

The  simplest  conceptual  and  physical  picture  of  the  evo¬ 
lution  from  BCS  to  BEC  superfluidity  for  s-wave  pairing  can 
be  constructed  for  low  fermion  densities  and  short-ranged 
interactions,  for  which  the  interaction  range  is  much  smaller 
than  the  average  separation  between  fermions.  In  that  limit 
the  essence  of  the  evolution  at  zero  temperature  is  reflected 
in  the  ratio  of  the  size  of  fermion  pairs  to  the  average  sepa¬ 
ration  between  the  fermions.  In  the  BCS  regime,  the  attrac¬ 


tion  is  weak,  the  pairs  are  much  larger 
than  their  average  separation,  and  they 
overlap  substantially.  In  the  BEC  regime, 
the  attraction  is  strong,  the  pairs  are 
much  smaller  than  their  average  separa¬ 
tion,  and  they  overlap  very  little. 

Amazingly,  a  clear  picture  of  the 
BCS-to-BEC  evolution  at  zero  tempera¬ 
ture  didn't  emerge  until  1980,  when  An¬ 
thony  Leggett  realized  that  the  physics 
could  be  captured  by  a  simple  descrip¬ 
tion  in  real  space  of  paired  fermions  with 
opposite  spins.5  Leggett  considered  a 
zero-ranged  attractive  potential  — that  is, 
a  contact  interaction— between  fermions 
and  showed  that  when  the  attraction  was 
weak,  a  BCS  superfluid  appeared,  and 
when  the  attraction  was  strong,  a  BEC  su¬ 
perfluid  emerged.  Philippe  Nozieres  and 
Stephan  Schmitt-Rink  (NSR)  used  a  dia¬ 
grammatic  method  with  a  finite-ranged 
attraction  to  extend  Leggett's  description 
to  temperatures  near  the  critical  temper¬ 
ature  for  superfluidity.6  Such  ideas  re¬ 
mained  largely  academic,  however,  since 
the  evolution  from  the  BCS  to  the  BEC 
regime  had  not  materialized  in  the  ex¬ 
perimental  world. 

But  the  discovery  of  cuprate  high-T. 
superconductors  in  1986  created  a  new 
paradigm  for  superconductivity.  The 
BCS  theory  seemed  to  fail  dramatically  in 
some  important  regions  of  the  phase  di¬ 
agram  of  cuprate  superconductors,  some 
of  which  have  critical  temperatures  near 
100  K.  Motivated  by  the  inapplicability  of 
the  BCS  theory  to  cuprate  superconduc¬ 
tors,  which  seemed  to  have  small  electron  pairs,  Jan  Engel- 
brecht,  Mohit  Randeria,  and  I  extended  the  preliminary  re¬ 
sults  of  Leggett  and  NSR  as  an  attempt  to  understand 
cuprates.7  We  used  a  zero-ranged  attraction  characterized  by 
the  experimentally  measurable  length  scale  the  so-called 
scattering  length.  The  natural  momentum  scale  is  the  Fermi 
wavenumber  kr,  and  the  strength  of  the  attractive  interactions 
can  be  described  by  the  dimensionless  scattering  parameter 
1  /kfa„.  That  parameter  changes  sign  in  going  from  weak  to 
strong  attraction:  The  BCS  weak-attraction  limit  is  character¬ 
ized  by  l/kra<  «:  -1;  the  BEC  strong-attraction  limit,  by 
l/fcFfls  3>  +1.  The  crossover  region  between  the  two  limits  oc¬ 
curs  for  -1  <  Vkfa,  <+l. 

The  phase  diagram  in  figure  2  illustrates  two  important 
physical  concepts  of  the  evolution  from  BCS  to  BEC  super¬ 
fluidity.  First,  the  normal  state  for  weak  attractions  is  a  Fermi 
liquid,  which  evolves  smoothly  into  a  molecular  Bose  liquid. 
The  two  regimes  are  separated  by  a  pair  formation  (or  mo¬ 
lecular  dissociation)  temperature  T„lr  characterized  by  chem¬ 
ical  equilibrium  between  bound  fermion  pairs  and  unbound 
fermions.  Second,  T  ir  and  the  critical  temperature  Tc  are  es¬ 
sentially  the  same  in  the  BCS  limit:  Pairs  form  and  develop 
phase  coherence  (condense)  at  the  same  temperature.  But  in 
the  BEC  limit,  fermion  pairs  (diatomic  molecules)  form  first 
around  T  ir  and  condense  at  the  much  lower  temperature 
Tc  =  Tbec,  where  phase  coherence  is  established. 

Two  other  concepts  are  shown  in  figure  2.  First  is  the  so- 
called  unitarity  limit,  where  the  scattering  length  diverges. 
As  the  scattering  parameter  l/fcF<rs  changes  sign  from  nega¬ 
tive  to  positive,  a  two-body  bound  state  with  characteristic 
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Figure  3.  Evidence  of  superfluidity,  (a)  Two  maps  of  the  density  of  potassium-40 
in  the  regime  where  the  atoms  form  small  fermion  pairs,  or  molecular  bosons.  On 
the  left,  the  temperature  is  above  the  critical  temperature  and  the  molecules  are 
uncondensed.  On  the  right,  the  temperature  is  below  the  critical  temperature  and 
the  molecules  undergo  Bose-Einstein  condensation,  as  evidenced  by  the  sharp 

Ek.  (Adapted  from  ref.  10.)  (b)  More  striking  evidence  for  superfluidity,  here  of 
um-6.  Absorption  images  taken  in  the  crossover  stable  region  between  the 
Bardeen-Cooper-Schrieffer  and  BEC  regimes  (see  figure  2)  reveal  the  existence  of 
a  lattice  of  vortices.  (Adapted  from  ref.  1 1 .) 


size  equal  to  the  scattering  length  emerges  in  vacuum  (see  the 
box  on  page  47).  Yet  despite  the  vanishing  of  the  scattering 
parameter  and  the  emergence  of  a  molecular  bound  state, 
nothing  really  dramatic  happens  to  the  superfluid  state  at 
unitarity.  What  is  more  relevant  from  a  many-body  perspec¬ 
tive  is  the  vanishing  of  the  chemical  potential  (/x  =  0),  which 
occurs  beyond  the  unitarity  limit  at  zero  temperature.  That 
special  point  is  the  de  facto  separation  between  the  BCS  re¬ 
gion,  where  the  energy  gap  in  the  elementary  excitation  spec¬ 
trum  is  related  to  the  order  parameter  and  occurs  at  finite  mo¬ 
mentum,  and  the  BEC  region,  where  the  energy  gap  occurs 
at  zero  momentum  and  is  related  to  both  the  chemical  po¬ 
tential  and  the  order  parameter.  It  is  the  qualitative  change 
in  the  elementary  excitation  spectrum  at  zero  chemical  po¬ 
tential— not  the  emergence  of  a  bound  state  in  vacuum  at  the 
unitarity  limit,  where  the  scattering  length  diverges— that 
separates  the  BCS  region  from  the  BEC  region  at  zero  tem¬ 
perature.  But  that  raises  the  question.  Is  it  possible  to  tune  the 
interactions  to  cross  the  line  of  zero  chemical  potential  and 
go  from  one  region  to  the  other? 

Tuning  the  interactions 

In  superconductors  some  tunability  of  the  fermion  density 
can  be  achieved  through  chemical  doping  or  electrostatic  gat¬ 
ing  in  the  same  material,  and  in  liquid  ’He  some  change  in 
density  is  possible  under  pressure,  but  essentially  no  control 
over  interactions  is  possible.  In  nuclei  and  neutron  stars,  the 
situation  is  even  worse— there  is  no  control  at  all  over  the 
density  or  interaction  strength.  But  five  years  ago,  the  tuning 
of  interactions  became  possible  in  ultracold  Fermi  atoms 
through  magnetically  driven  Feshbach  resonances  (see  the  box 
and  also  Dan  Kleppner's  column,  PHYSICS  TODAY,  August  2004, 
page  12).  It  is  this  ability  that  makes  ultracold  atoms  a  special 
laboratory  to  study  strongly  correlated  fermion  systems. 

In  late  2003,  almost  simultaneously,  three  experimental 
groups  succeeded  in  producing  BEC  of  bound  fermion  mol¬ 
ecules  in  optical  traps  at  ultracold  temperatures.  Two  groups, 
Rudolf  Grimm's  at  the  University  of  Innsbruck8  and  Wolf¬ 
gang  Ketterle's  at  MIT,9  used  6Li  atoms,  and  a  third  group, 
Deborah  Jin's  at  JILA,10  used  40K.  Absorption  images  of  the 
atoms  after  they  were  released  from  the  trap  showed  a  sharp 
central  peak  in  the  density  (figure  3a),  an  important  signa¬ 
ture  of  the  condensation  of  molecular  bosons  (that  is,  tightly 
paired  fermions).  But  the  true  evidence  for  superfluidity  was 
the  experimental  detection  by  the  MIT  group"  of  a  vortex  lat¬ 


tice  in  6Li  when  the  cloud  of  trapped 
fermions  was  first  rotated  and  then 
allowed  to  expand  (see  figure  3b). 
Vortices  are  stable  topological  excita¬ 
tions  that  are  formed  due  to  circulat¬ 
ing  superfluid  currents  and  carry  one 
unit  of  angular  momentum.  The  vor¬ 
tices  interact  repulsively  to  form  a  tri¬ 
angular  lattice  that  is  preserved  even 
upon  expansion  of  the  cloud  since  the 
topological  excitations  are  quite  stable. 

Additional  experimental  devel¬ 
opments  at  MIT12  and  Rice  Univer¬ 
sity13  in  2005  initiated  the  study  of 
Fermi  systems  (6Li)  with  population 
imbalance— that  is,  having  unequal 
numbers  of  fermions  in  distinct  hy- 
perfine  states,  labeled  "up"  and 
"down."  Starting  from  an  equal  pop¬ 
ulation  mixture,  appropriate  applica¬ 
tion  of  RF  pulses  can  arbitrarily 
change  the  relative  populations  by 
converting  up  fermions  to  down  and  vice  versa.  The  ability 
to  control  the  population  of  fermions  and  create  imbalances 
varying  almost  continuously  from  equal  mixtures  to  only  one 
occupied  state  has  stimulated  substantial  experimental  and 
theoretical  work  to  understand  phase  diagrams  throughout 
the  evolution  from  BCS  to  BEC  superfluidity  at  zero  and  fi¬ 
nite  temperatures. 

One  expects  a  minimum  of  three  phases  when  the  pop¬ 
ulation  imbalance  is  fixed  and  the  interaction  is  changed.  As 
the  interaction  increases  at  zero  temperature,  the  system  is 
first  a  normal,  unpaired  fluid.  It  then  separates  into  distinct 
regions,  one  containing  a  paired  superfluid  and  the  other 
containing  excess  unpaired  fermions,  before  it  finally  reaches 
a  phase  in  which  the  superfluid  and  excess  fermions  coexist 
(see  figure  4).  A  key  feature  of  the  phases  is  that  they  are  per¬ 
fectly  symmetric  with  respect  to  population  imbalance,  and 
it  does  not  matter  if  the  excess  fermions  are  up  or  down,  since 
atoms  have  the  same  mass. 

The  next  frontiers 

The  unprecedented  control  and  tunability  in  ultracold  fermi¬ 
ons  may  illuminate  the  physics  of  strongly  correlated  fermi¬ 
ons  in  standard  condensed-matter  physics,  nuclear  physics, 
and  astrophysics,  where  control  is  much  more  limited  or  non¬ 
existent.  Among  the  many  possible  experimental  and  theoret¬ 
ical  directions,  five  are  perhaps  of  the  most  immediate  nature. 
►  Mixtures  of  fermions  of  unequal  masses.  New  possibili¬ 
ties  arise  when  one  considers  mixtures  of  fermions  of  un¬ 
equal  masses,  such  as  6Li  and  40K,  '’Li  and  s7Sr,  or 40 K  and  87Sr. 
Studying  such  mixtures  is  significant  not  only  for  atomic 
physics,  but  also  for  condensed-matter  physics,  where  they 
could  simulate  fermions  with  different  effective  masses,  and 
for  quantum  chromodynamics  (QCD),  where  they  could  sim¬ 
ulate  mixtures  of  quarks  with  unequal  masses.  As  in  the  case 
of  equal  masses,  a  minimum  of  three  phases  (normal,  phase- 
separated,  and  coexisting)  is  expected  for  unequal  masses 
throughout  the  evolution  from  BCS  to  BEC  superfluidity. 
However,  there  is  a  dramatic  asymmetry  in  the  phase  dia¬ 
gram  of  population  imbalance  versus  interaction  parameter. 

Conceptually,  the  easiest  limit  to  understand  is  the  BEC 
regime  of  a  mixture  of  equal  numbers  of  up  and  down  fermi¬ 
ons.  If  their  masses  are  different,  they  form  heteronuclear  di¬ 
atomic  molecules  that  repel  each  other  weakly  and  produce 
a  molecular  BEC.  When  excess  fermions  are  present,  the 
paired  fermions  also  repel  the  excess  fermions,  but  the  re- 


48  October  2008  Physics  Today 


www.physicstoday.org 


Normal  Phase  separation  Coexistence 


Figure  4.  Population  imbalances,  which  can  be  created  by  applying  RF 
pulses,  can  affect  the  superfluid  state  as  it  evolves  from  the  Bardeen-Cooper- 
Schrieffer  to  the  Bose-Einstein  condensation  regime.  Some  of  the  possible 
phases  at  zero  temperature  are  illustrated  here  for  an  excess  of  "up"  (blue) 
atoms.  For  a  fixed  and  sufficiently  large  population  imbalance  the  ground 
state  of  the  system  with  weak  attraction  between  the  fermions  is  a  normal 
mixture  of  unpaired  "up"  and  "down"  (red)  atoms— they  cannot  overcome  a 
pairing  barrier  created  by  their  population  imbalance.  As  the  attraction  in¬ 
creases,  the  blue  and  red  atoms  acquire  enough  energy  to  pair  with  each 
other  but  repel  the  excess  atoms,  which  leads  to  phase  separation  between 
the  superfluid's  paired  fermions  and  excess,  unpaired  fermions.  As  the  at¬ 
traction  between  fermions  is  increased  further,  the  pairs  become  more  tightly 
bound  and  repel  the  excess  atoms  less;  the  result  is  a  mixed  state  of  coexist¬ 
ing  superfluid  and  excess  fermions. 


pulsion  is  smaller  when  lighter  fermions  are 
in  excess,  and  larger  when  heavier  fermions 
are  in  excess,  which  makes  the  coexistence 
region  between  pairs  and  lighter  excess 
fermions  larger  than  that  of  pairs  and  heav¬ 
ier  excess  fermions.  That  asymmetry  awaits 
experimental  confirmation,  but  two  major 
experimental  advances  are  narrowing  the 
gap  with  theory:  Quantum  degenerate  mix¬ 
tures  of  6Li  and  40 K  have  been  created  ex¬ 
perimentally,14  and  Feshbach  resonances 
that  permit  the  control  of  interactions  in 
such  mixtures  have  been  reported.15 

►  BCS-to-BEC  evolution  for  higher  angu¬ 
lar-momentum  pairing.  Although  the  evo¬ 
lution  from  BCS  to  BEC  superfluidity  for 
s-wave  pairing  is  fairly  well  understood 
theoretically  and  experimentally,  there  is  no 
experimental  realization  of  the  phenome¬ 
non  for  p- wave  pairing.  Even  though 
p- wave  pairing  has  been  extensively  stud¬ 
ied  in  connection  with  superfluid  liquid 
3He,  some  superconductors,  and  even  neu¬ 
tron  stars,  most  theoretical  work  has  built 
on  the  BCS  theory,  under  the  assumption 
that  the  systems  are  somewhat  weakly  in¬ 
teracting.  However,  ultracold  fermions  like 
6Li  and  40K  have  some  known  p- wave  Fesh¬ 
bach  resonances  between  atoms  in  identical 
hyperfine  states,  and  those  may  allow  the 
tuning  of  p- wave  interactions. 

Although  p- wave  resonances  are  typi¬ 
cally  much  narrower  than  s-wave  and  are  much  more  diffi¬ 
cult  to  study,  they  possess  much  richer  physics  due  to  their 
anisotropy.  A  simple  conceptual  example  of  that  richness  is 
the  presence  of  an  anisotropic  elementary  excitation  spec¬ 
trum  that  is  gapless  in  the  BCS  regime  but  becomes  fully 
gapped  in  the  BEC  regime.  That  evolution  represents  a  phase 
change  that  is  not  accompanied  by  a  change  in  the  order- 
parameter  symmetry  (p-wave  in  both  regimes),  and  thus  it 
cannot  be  fully  described  by  current  theories  of  phase  tran¬ 
sitions.  An  experimental  observation  of  the  phenomenon 
would  be  a  textbook  example  of  the  need  for  a  new,  general¬ 
ized  theory. 

►  Effects  of  disorder  in  BCS-to-BEC  evolution.  Although 
disorder  is  difficult  to  control  in  standard  condensed-matter 
systems,  random  potentials  for  trapped  ultracold  fermions 
can  be  produced  under  controlled  circumstances.  That  allows 
the  study  of  the  effects  of  weak  and  strong  disorder  through¬ 
out  the  evolution  from  BCS  to  BEC  superfluidity  and  the 
study  of  three-dimensional  phase  diagrams  involving  tem¬ 
perature,  interaction,  and  disorder. 

To  give  a  sense  of  how  rich  that  phase  diagram  can  be, 
consider  first  a  random  potential  that  is  independent  of  the 
hyperfine  state  of  the  atoms.  Such  a  potential  corresponds  to 
the  case  of  nonmagnetic  impurities  in  standard  condensed- 
matter  systems.  For  an  s-wave  BCS  superfluid,  the  amplitude 
and  phase  of  the  order  parameter  are  strongly  coupled,  such 
that  the  breaking  of  Cooper  pairs  and  destruction  of  phase 
coherence  occur  simultaneously.  Because  the  random  poten¬ 
tial  is  not  pair-breaking,  it  does  not  affect  the  phase  coher¬ 
ence  associated  with  the  order  parameter,  and  thus  the  criti¬ 
cal  temperature  in  the  presence  of  weak  disorder  is 
essentially  unchanged.  The  robustness  of  s-wave  BCS  super¬ 
fluids  with  time-reversed  fermion  pairing  (for  example,  pairs 
with  opposite  momenta  and  spin)  in  the  presence  of  weak 
disorder  was  proposed  nearly  50  years  ago  by  Philip  Ander¬ 


son  and  is  known  as  Anderson's  theorem.  For  weakly  disor¬ 
dered  BEC  superfluids,  in  contrast,  the  breaking  of  local  pairs 
and  the  loss  of  phase  coherence  occur  at  very  different  en¬ 
ergy  scales.  The  critical  temperature  is  strongly  affected  by 
weak  disorder,  since  phase  coherence  is  more  easily  de¬ 
stroyed  without  the  need  to  break  local  pairs  simultaneously, 
and  Anderson's  theorem  doesn't  apply. 

One  can  also  study  the  phase  diagram  for  s-wave  super¬ 
fluids  in  random  potentials  that  are  dependent  on  the  hyper¬ 
fine  state  of  the  atoms;  such  potentials  simulate  the  effects  of 
magnetic  impurities  in  standard  condensed  matter.  The  effects 
of  weak  disorder  in  the  BCS  regime  in  this  case  are  substan¬ 
tial,  since  the  random  potential  is  pair-breaking,  which  leads 
to  the  loss  of  phase  coherence  and  to  a  quick  suppression  of 
Tc.  The  Tc  reduction  in  the  BEC  regime  is  also  strong,  as  with 
state-independent  random  potentials,  since  phase  coherence  is 
again  quickly  destroyed  without  pair  breaking.  However,  the 
effects  of  strong  disorder  on  superfluidity  are  more  complex, 
thanks  to  disorder-induced  localization— termed  Anderson  lo¬ 
calization— of  fermions  in  the  normal  state  of  the  BCS  regime 
and  of  bosons  in  the  normal  state  of  the  BEC  regime.  The  un¬ 
derstanding  of  the  normal  and  superfluid  states  in  the  strong- 
disorder  limit  throughout  the  evolution  from  BCS  to  BEC  is 
still  a  theoretical  challenge,  but  experimentalists  may  be  able 
to  reach  that  exciting  regime  in  a  controlled  way  and  to  study 
disorder-induced  phase  transitions. 

►  Mixtures  of  fermions  with  three  hyperfine  states.  Ex¬ 
panding  the  number  of  hyperfine  states  from  two  to  three, 
which  can  be  labeled  red,  blue,  and  green  (or  1, 2,  and  3),  opens 
quite  interesting  possibilities,  as  there  can  be  three  types  of  s- 
wave  pairings:  red-blue,  blue-green,  and  green-red.  Such 
pairs  may  form  a  color  superfluid  with  the  coexistence  of  the 
three  types  of  superfluids.  That  is  analogous  to  a  situation 
encountered  in  QCD  in  which  quarks  may  pair  in  different 
color  states  in  the  cores  of  neutron  stars  and  produce  a  color 
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Figure  5.  Simulating  other  systems.  The  ability  to  control  interactions  among  ultracold  fermions  allows  connections  to  be 
made  with  other  areas  of  physics  and  may  help  address  some  of  the  key  issues  involving  strongly  correlated  fermions  in 
neutron  stars,  nuclei,  and  condensed  matter,  (a)  Quarks  of  color  quantum  states  red,  blue,  and  green  may  pair  at  "low" 
temperatures  and  high  densities,  producing  a  color  superconductor  in  the  core  of  a  neutron  star.  Color  superconductivity 
might  be  simulated  in  ultracold  Fermi  atoms  if  a  stable  mixture  of  three  hyperfine  states  for  the  same  type  of  atom  can  be 
produced,  (b)  Cuprate  superconductors,  which  have  this  complicated  phase  diagram,  might  be  simulated  by  ultracold 
fermions  in  optical  lattices.  Although  the  temperature  scales  and  the  densities  for  ultracold  atoms  are  much  smaller  than  the 
corresponding  temperatures  and  densities  for  color  or  cuprate  superconductors,  the  ratios  of  the  temperature  to  the  Fermi 
temperature  may  be  close,  permitting  some  parallels  to  be  drawn.  It  would  be  remarkable  to  achieve  color  superconductiv¬ 
ity  without  quarks,  and  high-temperature  superconductivity  in  a  lattice  without  copper  and  oxygen. 


superconductor  at  densities  several  times  larger  than  those 
found  in  typical  nuclear  matter  (see  figure  5a). 

Important  differences  exist,  however,  between  the  color 
superconductor  expected  in  QCD  and  the  color  superfluid 
expected  in  ultracold  fermions.  First,  ultracold  atoms  are 
neutral  and  quarks  are  fractionally  charged  particles.  Second, 
the  densities  of  ultracold  fermions  are  extremely  low,  while 
the  densities  for  quarks  are  extremely  high.  Third,  whereas 
quarks  have  a  threefold  degeneracy  of  their  color,  the  de¬ 
generacy  of  hyperfine  states  in  ultracold  atoms  is  usually 
lifted  by  an  external  magnetic  field.  The  need  to  control  two- 
body  collisions  and  prevent  three-body  collisions  of  trapped 
atoms  involving  three  distinct  hyperfine  states  makes  exper¬ 
iments  difficult,  and  color  superfluidity  has  not  yet  material¬ 
ized  in  ultracold  fermions.  Still,  it  is  not  crazy  to  think  about 
it!  A  natural  candidate  for  color  superfluidity  is  6Li,  which  has 
three  hyperfine  states  that  might  be  suitable,  if  losses  can  be 
cleverly  controlled.  Other  ultracold  fermions,  including 
Fermi  isotopes  of  ytterbium,  are  on  the  horizon.  Furthermore, 
there  is  the  promise  of  being  able  to  control  the  interactions 
and  populations  of  ultracold  atoms  to  obtain  and  explore 
phases  that  do  not  yet  have  analogous  behavior  in  QCD. 

►  Ultracold  fermions  in  optical  lattices.  The  evolution  from 
BCS  to  BEC  superfluidity  for  s-wave  systems  in  3D  cubic  lat¬ 
tices  was  first  discussed  by  NSR  in  1985,  but  experiments 
probing  superfluidity  of  ultracold  fermions  in  optical  lattices 
are  just  beginning.  Tilman  Esslinger  and  colleagues  at  ETH 
Zurich  reported  the  production  of  “K  molecules  in  3D  cubic 
optical  lattices  using  s-wave  Feshbach  resonances  in  early 
2006,16  but  no  evidence  of  a  superfluid  state  was  found  until 
later  that  year,  when  Ketterle  and  coworkers  loaded  6Li  atoms 
in  optical  lattices  and  their  pairs  formed  a  condensate.17 
Those  two  experiments  opened  the  door  to  studies  of  super¬ 


fluid-to-insulator  transitions  in  optical  lattices  and  are  likely 
to  stimulate  and  renew  the  interest  on  the  possible  realiza¬ 
tion  of  d- wave  superfluidity  in  nearly  2D  lattices  like  those 
encountered  in  the  high-Tc  cuprate  superconductors,  which 
have  a  complex  phase  diagram  (see  figure  5b)  and  are  known 
to  be  d- wave  superconductors.  It  would  be  quite  interesting 
to  tune  interactions  and  filling  factors  of  ultracold  fermions 
in  optical  lattices,  and  try  to  reproduce  the  phase  diagram  of 
cuprate  superconductors. 

Predicting  the  future 

The  unprecedented  control  over  interactions  in  various  di¬ 
mensionalities  and  geometries  has  put  research  of  ultracold 
fermions  in  harmonic  traps  and  optical  lattices  at  the  fore¬ 
front  of  investigations  into  the  behavior  of  Fermi  condensates 
and  strongly  interacting  fermions.  In  particular,  exotic  phases 
of  QCD,  like  color  superconductivity,  and  superfluid  phases 
in  neutron  stars  and  nuclear  matter  may  have  analogous 
counterparts  in  tabletop  experiments  involving  ultracold 
fermions.  Furthermore,  studies  of  disorder-induced  phe¬ 
nomena  are  also  within  experimental  reach,  and  insights  into 
p-  and  d- wave  superconductivity  found  in  condensed-matter 
physics  are  just  around  the  corner,  as  ultracold  fermions 
loaded  into  optical  lattices  begin  to  be  explored.  No  one  has 
a  crystal  ball  to  predict  the  future,  but  research  on  ultracold 
atoms  will  likely  provide  insight  into  the  simulation  and  un¬ 
derstanding  of  many  known  and  unknown  phases  of  nuclear, 
atomic,  molecular,  and  condensed-matter  physics. 

I  thank  Sam  Bader,  Ian  Spielntan,  Kris  Helmerson,  and  Steve  Anlage, 
among  other  colleagues,  for  reading  this  manuscript  at  various  stages 
of  its  production  and  for  their  valuable  suggestions.  This  article  is  ded¬ 
icated  to  my  father  and  mother,  who  passed  away  during  its  writing. 
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Preparing  a  highly  degenerate  Fermi  gas  in  an  optical  lattice 
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(Dated:  November  6,  2008) 

We  propose  a  method  to  prepare  a  sample  of  fermionic  atoms  in  a  three-dimensional  (3D)  optical 
lattice  at  unprecedentedly  low  temperatures  and  uniform  filling  factors.  The  process  involves  adia¬ 
batic  loading  of  atoms  into  multiple  energy  bands  of  an  optical  lattice  followed  by  a  filtering  stage 
whereby  atoms  from  all  but  the  ground  band  are  removed.  Of  critical  importance  is  the  use  of  a  non¬ 
harmonic  trapping  potential,  taken  here  to  be  the  radial  profile  of  a  high-order  Laguerre-Gaussian 
laser  beam,  to  provide  external  confinement  for  the  atoms.  For  realistic  experimental  parameters, 
this  procedure  should  produce  samples  with  temperatures  ~  10-3  of  the  Fermi  temperature.  This 
would  allow  the  investigation  of  the  low-temperature  phase  diagram  of  the  Fermi-Hubbard  model 
as  well  as  the  initialization  of  a  high-fidelity  quantum  register. 

PACS  numbers:  03.75.Ss,  32.80.Pj,  03.67.Lx,  37.10.De,  37.10.Jk,  05.30.-d 


Investigations  of  degenerate  Fermi  gasses  loaded  into 
optical  lattices  have  indicated  that  these  systems  are 
ideal  for  creating  a  robust  quantum  register  for  quan¬ 
tum  computing  applications  [1,  2]  as  well  as  providing 
a  testing  ground  for  paradigm  models  of  condensed  mat¬ 
ter  physics.  Models  currently  under  investigation  include 
studying  Fermi  surfaces  and  band  insulator  states  [3], 
fermionic  superfluidity  in  a  lattice  [4],  and  transport 
properties  of  interacting  fermions  in  one  and  three  dimen¬ 
sional  optical  lattices  [5,  6].  These  seminal  experiments 
demonstrate  the  high  precision  and  versatility  available 
in  simulating  solid  state  systems  with  fermions  in  optical 
lattices. 

Theoretical  studies  of  such  systems  have  predicted  that 
a  number  of  exotic  phases  emerge  at  low  temperatures, 
including  quantum  magnetic  ordering  and  possibly  d- 
wave  superfluidity  [7,  8].  However,  temperatures  low 
enough  to  observe  exotic  phases  such  as  these  are  dif¬ 
ficult  to  achieve  when  optical  lattices  are  loaded  with  a 
harmonic  external  confining  potential.  It  has  been  the¬ 
oretically  predicted  [9]  and  experimentally  observed  [3] 
that  fermions  adiabatically  loaded  into  an  optical  lattice 
with  harmonic  external  confinement  experience  heating 
for  all  but  very  high  initial  temperatures  and  filling  fac¬ 
tors  (number  of  atoms  per  lattice  site)  [10]. 

Alternative  methods  to  prepare  fermionic  atoms  in  op¬ 
tical  lattices  at  low  temperatures  and/or  high  uniform 
filling  factors  include:  cooling  by  adiabatic  loading  into 
a  three-dimensional  (3D)  homogeneous  trapping  poten¬ 
tial  with  high  filling  factor  [11],  defect  filtering  in  a  state 
dependent  optical  lattice  [1],  adiabatic  loading  [2]  and  fil¬ 
tering  [12]  of  high  entropy  atoms  from  a  ID  lattice  with 
harmonic  confinement. 

In  this  Letter,  we  propose  a  method  to  prepare  a  highly 
degenerate  Fermi  gas  in  a  3D  optical  lattice  using  a  box¬ 
like  potential  for  external  confinement  and  taking  advan¬ 
tage  of  the  Pauli  exclusion  principle  to  selectively  remove 
atoms  from  multiply-occupied  lattice  sites.  Specifically, 
we  assume  that  the  radial  profile  of  a  blue-detuned,  high- 


FIG.  1:  (color  online)  We  filter  high  entropy  atoms  from  a 
combined  box-like  potential  (blue  dashed  line)  and  optical 
lattice  (red  solid  line)  by  selectively  removing  atoms  from 
all  but  the  ground  energy  band.  Amplitude  modulation  of 
the  lattice  potential  can  selectively  transfer  these  atoms  to 
high-lying  bands  via  a  two-photon  transition  where  they  then 
tunnel  out  of  the  region.  Dramatic  cooling  results  when  the 
Fermi  energy  (prior  to  filtering)  lies  within  the  second  band. 


order  Laguerre-Gaussian  (LG)  laser  beam  provides  con¬ 
finement  along  each  cartesian  axis.  The  atoms  are  pre¬ 
pared  via  a  two  step  process:  (1)  adiabatically  loading 
atoms  initially  confined  in  the  LG  trap  into  a  superim¬ 
posed  optical  lattice,  followed  by  (2)  irreversibly  filtering 
atoms  from  all  but  the  ground  energy  band  (see  Fig.  1). 
We  find  that  when  the  Fermi  energy  of  the  system  is 
sufficiently  large,  such  that  atoms  begin  to  significantly 
populate  the  second  energy  band  prior  to  filtering,  con¬ 
siderable  cooling  is  achieved;  whereas  for  lower  filling  fac¬ 
tors  heating  is  observed. 

The  energy  spectrum  of  a  system  of  ultracold  atoms 
is  greatly  affected  by  the  addition  of  a  3D  cubic  optical 
lattice  which  can  be  formed  by  three  perpendicular  sets 
of  retroflected  Gaussian  laser  beams  detuned  far  from 
resonance.  In  a  homogeneous  trapping  potential,  the  lat¬ 
tice  breaks  the  translational  symmetry  of  the  system,  re¬ 
sulting  in  a  series  of  discrete  energy  bands  whose  width 
and  energy  spacings  are  dependent  on  the  intensity  and 
wavelength  of  the  laser  light  (see  Fig.  2(a)).  Bezett  and 
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FIG.  2:  (color  online)  (a)  Band  structure  for  a  5  Er  deep 
lattice,  (b)  Atoms  can  be  selectively  excited  from  the  n  to 
n  4-  2  energy  band  by  modulating  the  amplitude  of  the  op¬ 
tical  lattice  light.  The  excitation  energy  between  bands  1 
and  3  (blue)  and  2  and  4  (red)  for  fermions  in  all  crystal 
momenta  states  are  shown.  As  the  depth  of  the  optical  lat¬ 
tice  approaches  35  Er,  the  excitation  energy  bands  separate, 
thereby  demonstrating  the  feasibility  of  performing  the  band 
selective  excitations  required  in  our  filtering  method. 


Blakie  demonstrated  that  this  energy  band  structure  can 
be  exploited  to  dramatically  increase  the  degeneracy  of 
the  sample  for  a  homogeneous  system  [11].  For  a  dense 
atomic  gas,  with  a  filling  factor  greater  than  unity,  ap¬ 
plication  of  the  lattice  increases  the  Fermi  energy,  since 
it  lies  within  the  first  excited  band,  and  compresses  the 
Fermi  surface  resulting  in  a  dramatic  reduction  in  the 
degeneracy  temperature  T/Tp,  where  Tp  is  the  Fermi 
temperature. 

This  band  structure  also  permits  state-selective  opera¬ 
tions  to  manipulate  and  probe  the  energy  distribution  of 
the  sample.  One  such  method  involves  modulating  the 
depth  of  the  optical  lattice  to  selectively  excite  atoms 
from  the  ntori+2  energy  band  with  no  change  in  the 
crystal  momentum  q  [13].  In  contrast  to  non-interacting 
Bose  systems,  where  a  single  q  can  be  macroscopically  oc¬ 
cupied,  a  Fermi  system  necessarily  begins  to  fill  the  bot¬ 
tom  band  and  A q  =  0  transitions  must  be  excited  for  all 
occupied  values  of  q.  In  Fig.  2(b),  we  show  band  excita¬ 
tion  energies  as  a  function  of  lattice  depth  for  n  =  1  — >  3 
and  n  =  2  — *  4  transitions  spanning  all  q  within  a  Bril- 
louin  zone.  By  loading  the  sample  into  an  optical  lattice 
with  a  depth  of  Vo  =  35  Er  (where  Er  =  h2k2 /‘2m  is 
the  recoil  energy  and  k  is  the  wavenumber  of  the  laser 
light),  we  find  that  these  transitions  are  well  resolved. 
It  is  therefore  possible  to  apply  a  filtering  process  which 
selectively  removes  atoms  from  all  but  the  ground  energy 
band.  Using  adiabatic  rapid  passage,  population  may  be 
selectively  transferred  from  n  =  2  — *  4  by  sweeping  the 
amplitude  modulation  frequency  from  below  to  above  all 
2  — »  4  transition  frequencies  while  remaining  below  the 
lowest  1  — ♦  3  transition  frequency.  Then,  lowering  the 
height  of  the  trapping  potential  allows  atoms  in  the  third 
and  higher  energy  bands  to  tunnel  out  of  the  system. 

In  order  to  experimentally  approximate  the  homoge¬ 
nous  lattice  potential  described  above,  we  consider  the 


addition  of  a  box-like  external  potential  produced  by  a 
blue-detuned,  £‘h-order,  LG  laser  beam  with  a  radial  pro¬ 
file 


VLG(r)  -  Vpeak  (  w2(  )  «  2r  /w°  (1) 

at  the  beam  waist  wq.  For  a  given  charge  t,  the  peak 
value  Vpeak  of  the  potential  occurs  at  rmax  =  u>0  \fij2 
and  the  width  of  this  peak  decreases  with  decreasing  wo- 
Therefore,  for  a  given  trap  size  rmax,  the  LG  profile  more 
closely  approximates  a  box  potential  when  wq  is  reduced 
and  t  is  correspondingly  increased.  Trapping  of  ultracold 
gases  has  been  demonstrated  in  single  or  crossed  beam 
configurations  of  LG  beams  up  to  £  =  16  [14-16]. 

Along  a  given  cartesian  axis,  we  take  the  single  particle 
Hamiltonian  to  be 


H(x)  = 


- h 2  d2 
2m  dx2 


+  VLG(z) 


(2) 


1 


+Vo  cos  2(kx  +  </>x)  +  -murar, 


where  the  third  term  represents  a  lattice  potential  of 
depth  Vo  and  phase  offset  <f>x.  We  also  include  a  har¬ 
monic  term  that  arises  if  red-detuned  Gaussian  beams 
are  used  to  produce  the  lattice  potential;  in  this  case 
ui  on  \/Vo.  The  ID  eigenvalues  and  eigenfunctions  for  a 
given  depth  of  the  optical  lattice  are  calculated  by  nu¬ 
merically  diagonalizing  the  Hamiltonian  (Eq.  2)  using 
the  method  described  in  [17].  For  sufficiently  shallow 
lattice  depths,  the  low  energy  eigenstates  are  delocalized 
and  closely  approximate  Bloch  states  in  the  first  band  of 
a  homogeneous  system.  However,  higher  energy  states 
are  either  localized  at  the  edges  of  the  trap  (i.e.  near 
x  =  rmax)  or  delocalized  and  correspond  to  Bloch  states 
in  higher  bands.  While  a  band  structure  picture  is  not 
strictly  valid  for  this  inhomogeneous  system,  we  classify 
the  set  of  eigenfunctions  without  nodes  to  constitute  the 
first  band,  Ciband,?** 

We  extend  this  model  to  three  dimensions  by  assuming 
a  separable  Hamiltonian  H30  =  H(x)  +  H(y)  +  H(z). 
For  simplicity,  we  assume  equal  lattice  depths  in  each 
direction.  The  3D  spectrum  (£m)  for  a  given  depth  of  the 
optical  lattice  is  then  generated  by  calculating  all  possible 
combinations  of  the  sum  £m  =  ej  +  tj  -t-e*  for  all  values  of 
the  ID  eigenenergies  (ep)  in  each  spatial  direction.  The 
3D  energy  spectrum  for  energy  states  in  the  first  band 
of  the  optical  lattice  is  calculated  in  a  similar  manner 

(£lband,m  =  elband,t  T  £  1  band ,j  T  <:  lband.k ) ■ 

In  calculating  thermodynamic  quantities  during  the 
proposed  cooling  method,  we  assume  constant  thermal 
equilibrium  before  and  after  the  selective  removal  of 
atoms  from  high-lying  bands.  Equilibrium  is  main¬ 
tained  by  elastic  collisions  in  a  50/50  mixture  of  spin- 1/2 
fermions  and  changes  in  the  trapping  potential  are  adi¬ 
abatic  with  respect  to  the  rethermalization  time  scale. 
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However,  we  also  assume  that  the  interactions  are  weak 
enough  to  not  significantly  modify  the  single-particle  en¬ 
ergy  spectrum,  £m.  We  therefore  use  £m  when  calculat¬ 
ing  the  following  quantities: 

N  =  2V _ l- _ 

„  1  +  exp[(£m  -  n)/kBT] ' 

E  =  251  +  “pI(£»"-m)/W  (3) 

i  ~  2  E  '“I'  +  “Pl<£»  -  c)/‘*JH  +  -  jfr". 

where  T  is  the  temperature,  n  is  the  chemical  potential 
of  an  atom  in  either  spin  state,  N  is  the  total  number  of 
atoms,  E  is  the  total  energy  in  the  system,  and  S  is  the 
total  entropy.  The  degeneracy  temperature  is  given  by 
T/Tp  where  kBTp  =  £jv,  the  energy  of  the  Nth  eigen¬ 
state  of  the  multi-band,  3D  spectrum.  After  filtering, 
the  thermodynamic  quantities  in  Eq.  3  are  calculated  for 
atoms  only  in  the  first  band  using  £iband,m- 

Our  proposed  method  for  cooling  the  atoms  is  com¬ 
prised  of  (1)  an  adiabatic  increase  in  the  lattice  depth 
starting  from  zero,  (2)  a  non-adiabatic  selective  filtering 
of  atoms  and  (3)  an  optional  adiabatic  change  to  a  final 
lattice  depth.  To  calculate  changes  in  thermodynamic 
quantities  during  these  stages  we  use  the  following  pro¬ 
cedures.  For  adiabatic  changes  of  the  potential  we  (1) 
calculate  S  for  a  given  N  and  initial  temperature  T\  us¬ 
ing  the  energy  spectrum  for  the  initial  potential  and  (2) 
numerically  solve  for  Hf  and  the  final  temperature  Tf  in 
Eq.  3  using  the  spectrum  for  the  final  potential,  assum¬ 
ing  N  and  S  are  conserved.  In  contrast,  for  the  selective 
filtering  stage  we  (1)  start  from  a  thermalized  sample  of 
Ni  atoms  at  temperature  T\  for  a  given  spectrum,  (2)  cal¬ 
culate,  given  this  distribution,  the  energy  Ef  and  number 
Nf  for  atoms  restricted  to  the  first  band,  and  (3)  solve  for 
Pf  and  the  temperature  Tf  in  Eq.  3  using  the  multi-band 
energy  spectrum  assuming  the  sample  equilibrates  with 
total  energy  Ef  and  number  Nf. 

We  consider  a  50/50  spin  mixture  of  6Li  atoms  ini¬ 
tially  trapped  in  a  LG  trapping  potential  with  £  =  12, 
kpeak  =  35  Er  and  rmax  =  13.5  pm.  For  reasonable 
lattice  beam  properties  ( k  =  27t/1064  nm  and  a  waist 
of  200  pm)  we  find  u  =  2n  (586Hz)  for  the  final  lattice 
depth  Vo,f  =  35  EB.  The  final  degeneracy  temperatures 
after  adiabatic  loading  and  filtering,  along  with  the  fi¬ 
nal  atom  number  are  shown  in  Fig.  3  for  various  initial 
degeneracy  temperatures  and  sample  sizes.  In  each  case 
<!>x  =  4>y  —  <t>z  —  0.  This  data  shows  that  the  thermo¬ 
dynamic  properties  of  the  system  are  highly  dependent 
on  the  initial  filling  factor  and  can  be  separated  into  two 
distinct  regions  A  and  B.  The  vertical  dashed  line  which 
separates  the  regions  represents  the  number  of  atoms  at 
which  the  Fermi  energy  enters  the  second  band. 

In  region  A,  the  Fermi  energy  before  filtering  lies  be¬ 
low  the  second  energy  band.  For  very  low  filling  factors, 


FIG.  3:  (color  online)  As  a  function  of  initial  atom  number 
we  report  (a)  the  final  degeneracy  temperature  and  (b)  the 
final  atom  number  after  implementing  the  proposed  cooling 
and  filtering  procedure  for  various  initial  temperatures  be¬ 
tween  0.005  and  0.3 Tf-  The  vertical  dashed  line  represents 
the  number  of  atoms  for  which  the  Fermi  energy  enters  the 
second  band.  The  trap  and  lattice  parameters  are  as  described 
in  the  text. 


an  increase  in  T/Tp  is  observed  for  all  initial  tempera¬ 
tures.  This  increase  in  T /Tf  occurs  because  Tf  decreases 
more  than  T  as  the  lattice  depth  increases.  Additionally, 
Fig.  3(b)  demonstrates  that  we  are  not  significantly  fil¬ 
tering  atoms  for  low  initial  temperatures.  As  the  Fermi 
energy  approaches  the  second  energy  band,  we  see  a  dra¬ 
matic  decrease  in  the  final  T/Tp ■  In  this  regime,  where 
atoms  are  beginning  to  significantly  occupy  the  localized 
energy  states  at  the  edge  of  the  trapping  potential,  a 
dramatic  increase  of  the  Fermi  energy  is  observed. 

In  region  B,  the  density  is  such  that  the  Fermi  energy 
before  filtering  lies  within  the  second  band.  In  this  re¬ 
gion  the  adiabatic  increase  of  the  lattice  depth  results  in 
a  dramatic  increase  in  Tp,  a  substantial  reduction  in  the 
temperature  T  [11],  and  allows  for  a  significant  reduc¬ 
tion  in  entropy  during  the  filtering  stage.  We  find  that 
significant  cooling  is  achieved  for  initial  temperatures  in 
the  vicinity  of  T\  =  0.1  Tf.  Above  this  initial  tempera¬ 
ture  cooling  is  less  efficient.  Below  this  initial  tempera- 
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ture,  the  final  T/Tp  after  filtering  saturates.  As  can  be 
seen  in  Fig.  3(b),  N(  is  extremely  insensitive  to  fluctua¬ 
tions  in  N-,  for  low  initial  temperatures.  For  example,  at 
T|  =  0.05  Tf,  a  variation  of  ±10%  around  Ni  =  1.6  x  105 
yields  a  variation  of  only  +0.09%/  —  0.2%  in  Nf. 

The  cooling  efficiency  and  number  filtering  were  de¬ 
pendent  on  the  choice  of  phases  (j>x ,  <j>y  and  4>z  due  to 
the  sensitive  effect  these  phases  had  on  the  the  location 
of  localized  edge  state  eigenenergies  relative  to  the  Fermi 
energy.  To  study  this  effect,  we  modeled  the  system  al¬ 
lowing  the  phase  in  each  direction  to  be  independently 
selected  from  the  set  <t>a  =  (0,  ir/10,  We  consid¬ 

ered  samples  with  an  initial  temperature  T\  =  0.05  Tf 
and  an  initial  number  N-,  =  1.6  x  10s  atoms,  parameters 
within  the  saturated  regime  for  all  choices  of  phase  and 
close  to  optimal  for  cooling.  From  the  set  of  all  possible 
phase  combinations,  we  find  an  average  final  tempera¬ 
ture  Tf  =  0.0031  Tf  where  10%  of  the  ensemble  achieved 
temperatures  below  T\0  =  0.0023  Tp  and  90%  were  below 
Tqq  =  0.004 Tf.  From  this  same  set,  we  find  an  average 
final  number  N(  =  1.20  x  105,  with  Nio  =  1.18  x  105, 
and  Nqo  =  1.22  x  105.  The  filtering  process  further  re¬ 
sults  in  a  substantial  reduction  in  entropy.  The  initial 
entropy  per  atom  s*  =  0.28  kp  is  reduced  to  an  average 
final  value  of  sf  =  0.024 with  «io  =  0.014  and 
S90  =  0.033  fcB. 

It  is  in  general  possible  to  prepare  atoms  at  a  low 
T/Tp  in  a  shallow  lattice  potential,  if  so  desired,  by  adi- 
abatically  reducing  the  lattice  depth  after  the  filtering 
stage.  Continuing  the  example  from  above,  when  the 
lattice  depth  is  reduced  to  5  Ep  we  find  an  average  fi¬ 
nal  temperature  Tf  =  0.002  Tf,  Ti0  =  0.0013  Tf,  and 
T90  =  0.0028  Tf. 

We  now  consider  the  effects  of  the  charge  t  of  the  LG 
beam  for  samples  with  an  initial  T;  =  0.05  Tf,  phase 
4>x  =  4>y  —  4>z  =  0,  final  lattice  depth  of  35  Ep,  and  var¬ 
ious  initial  atom  numbers.  For  each  f-value,  the  waist 
of  the  LG  beam  is  adjusted  such  that  the  number  of 
states  below  the  second  energy  band  is  held  constant  at 
1.22  x  105.  As  shown  in  Fig.  4,  the  cooling  efficiency  of 
this  procedure  is  highly  dependent  on  the  charge.  Note 
that  for  £  =  1,  which  approximates  harmonic  exter¬ 
nal  confinement,  the  final  degeneracy  temperature  Tf/Tp 
never  drops  below  its  initial  value  of  TJTf  —  0.05.  For 
t  8,  the  minimum  degeneracy  temperature  saturates 
to  Tf/Tp  <  0.003.  For  higher  values  of  £,  the  extent  of 
the  saturation  regime  grows.  We  believe  that  this  satu¬ 
ration  is  caused  by  localized  atoms  at  the  edges  of  the 
LG  potential  rethermalizing  into  higher  energy  bands. 

In  this  Letter  we  proposed  a  method  for  preparing  a 
sample  of  highly  degenerate  fermions  by  adiabatic  load¬ 
ing  into  a  combined  optical  lattice  and  “box-like”  trap¬ 
ping  potential  followed  by  selective  removal  of  atoms 
from  all  but  the  ground  energy  band.  Numerical  calcu¬ 
lations  for  sample  sizes  ~  10s  predict  that  temperatures 
~  10-3  Tf  can  be  prepared  in  this  manner.  This  method 


FIG.  4:  (color  online)  The  effects  of  the  charge  £  of  the 
Laguerre  -  Gaussian  trapping  potential  on  the  efficiency  of 
our  proposed  cooling  and  filtering  method.  For  each  data 
set,  the  initial  temperature  T,  =  0.05  Tf  and  the  phases 
<j>x  =  <f>y  =  <t>z  =  0.  For  all  £  values,  the  number  of  atoms 
at  which  the  Fermi  energy  enters  the  second  band  (vertical 
dashed  line)  is  held  constant. 


is  robust  against  initial  number  and  temperature  fluctu¬ 
ations  for  a  sufficiently  cold  initial  sample  of  atoms  and 
yields  samples  with  little  variance  in  the  final  number. 
While  the  selective  removal  of  atoms  must  occur  in  a 
deep  lattice  (in  order  to  spectrally  resolve  the  band  exci¬ 
tations),  subsequent  reduction  of  the  lattice  depth,  if  de¬ 
sired,  yields  a  modest  amount  of  additional  cooling.  We 
expect  that  this  method  can  be  scaled  to  larger  samples 
for  which  still  lower  degeneracy  temperatures  would  be 
attained  due  to  the  diminished  role  localized  edge  states 
would  play.  Further,  the  “box-like”  trapping  potential 
offers  an  ideal  spatial  profile  for  simulating  solid  state 
physics  with  degenerate  atoms  in  optical  lattices  as  the 
relatively  flat  central  region  allows  for  a  large  number 
of  delocalized  states  while  the  curvature  at  the  edges  of 
the  traps  removes  the  constraint  of  loading  exact  atom 
numbers  to  realize  insulating  states.  Atoms  prepared  in 
this  manner  should  be  sufficiently  cold  to  explore  quan¬ 
tum  spin  phases  of  fermionic  atoms  which  are  currently 
inaccessible,  and  could  provide  a  physical  realization  of 
an  essentially  perfect  quantum  register. 
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Effect  of  resonant  interactions  on  the  stability  of  a  three-state  Fermi  gas 
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We  investigate  the  stability  of  a  three  spin  state  mixture  of  ultracold  fermionic  ®Li  atoms  over 
a  range  of  magnetic  fields  encompassing  three  Feshbach  resonances.  For  most  field  values,  we 
attribute  decay  of  the  atomic  population  to  three-body  processes  involving  one  atom  from  each 
spin  state  and  find  that  the  three-body  loss  coefficient  varies  by  over  four  orders  of  magnitude.  We 
observe  high  stability  when  at  least  two  of  the  three  scattering  lengths  are  small,  rapid  loss  near  the 
Feshbach  resonances,  and  two  unexpected  resonant  loss  features.  At  our  highest  fields,  where  all 
pairwise  scattering  lengths  are  approaching  at  —  — 2140oo,  we  measure  a  three-body  loss  coefficient 
Ls  ~  5  x  10-22  cm6/s  and  a  trend  toward  lower  decay  rates  for  higher  fields  indicating  that  future 
studies  of  color  superfluidity  and  trion  formation  in  a  SU(3)  symmetric  Fermi  gas  may  be  feasible. 

PACS  numbers:  67.85.Lm,  34.50.-s,  67.85.-d,  03.75.Ss 


Multi-component  Fermi  gases  with  tunable  interac¬ 
tions  are  exceptionally  well  suited  to  the  study  of 
few-  and  many-body  quantum  physics.  Ultracold  two- 
state  Fermi  gases  near  a  Feshbach  resonance  have  been 
used  to  characterize  the  crossover  from  Bardeen-Cooper- 
Schreiffer  (BCS)  superfluidity  to  Bose-Einstein  condensa¬ 
tion  (BEC)  of  diatomic  molecules  [1-7].  Further,  a  nor¬ 
mal  to  superfluid  transition  and  phase  separation  have 
been  observed  in  imbalanced  two-state  spin  mixtures  [8- 
12].  The  stability  of  the  two-state  mixtures  against  two- 
and  three-body  loss  processes  was  critically  important  to 
the  success  of  these  experiments. 

Recently,  there  has  been  considerable  interest  in  the 
study  of  three-state  Fermi  gases  with  tunable  interac¬ 
tions  [13-25].  Whereas  three-body  interactions  in  ultra¬ 
cold  two-state  Fermi  gases  are  suppressed  by  the  exclu¬ 
sion  principle,  the  addition  of  a  third  spin  component 
allows  for  the  study  of  three-body  phenomena  such  as 
the  Efimov  effect  in  a  Fermi  system  if  pairwise  interac¬ 
tions  are  resonantly  enhanced  [19,  24].  Further,  with  this 
additional  spin  state  there  may  be  competition  between 
multiple  pairing  states  and  trion  formation  [13-17,  20- 
23].  If  pairwise  interactions  are  all  attractive  and  of  equal 
magnitude,  the  system  is  symmetric  under  SU(3)  trans¬ 
formations  and  is  predicted  to  exhibit  a  novel  superfluid 
phase  with  SU(2)  symmetry  and  an  anomolous  number 
of  Goldstone  modes  [15].  This  superfluid  phase  would  be 
a  cold  atom  analog  of  the  color  superfluid  phase  predicted 
in  quantum  chromodynamics  (QCD)  which  undergoes  a 
quantum  phase  transition  to  a  Fermi  liquid  of  trion  states 
(analogous  to  baryons)  as  the  ratio  of  the  interaction  en¬ 
ergy  to  the  kinetic  energy  is  adjusted  [18]. 

Future  studies  of  the  above  phenomena  critically  de¬ 
pend  on  the  magnitude  of  two-  and  three-body  loss  rates, 
particularly  when  two  or  more  scattering  lengths  are  en¬ 
hanced  by  scattering  resonances.  Two-  and  three-body 
loss  and  heating  processes,  often  enhanced  near  the  reso¬ 
nances  [26],  can  impose  stringent  limits  on  the  maximum 
achievable  phase  space  density.  Further,  the  unambigu¬ 


ous  observation  of  three-body  resonances  requires  negli¬ 
gible  two-body  loss  [27]. 

Most  investigations  of  Fermi  gases  with  tunable  inter¬ 
actions  have  focussed  on  two-state  mixtures.  A  small 
admixture  of  a  third  spin  component  has  been  used 
for  thermometry  [28]  and  recently  the  rapid  decay  of 
a  three-state  Fermi  gas  near  a  Feshbach  resonance  was 
noted  during  an  investigation  of  Cooper  pair  size  by 
radio-frequency  (RF)  spectroscopy  in  a  two-state  Fermi 
gas  [29].  Very  recently  our  group  reported  measurements 
of  three-body  loss  in  thermal  and  degenerate  three-state 
Fermi  gases  over  a  wide  range  of  magnetic  fields  [30]  and 
Ottenstein  et  al.  independently  reported  measurements 
of  three-body  loss  rate  coefficients  in  a  degenerate  three- 
state  Fermi  gas  for  field  values  away  from  any  pairwise 
scattering  resonances  [31]. 

In  this  Letter,  we  report  the  first  measurements  of 
three-body  loss  coefficients  for  a  three-state  Fermi  gas 
over  a  range  of  magnetic  fields  where  pairwise  interac¬ 
tions  are  resonantly  enhanced.  This  enhancement  is  due 
to  three  overlapping  interspecies  Feshbach  resonances 
and  a  zero-energy  resonance  in  the  6  Li  triplet  molecu¬ 
lar  potential.  Further,  we  demonstrate  that  all  two-state 
mixtures  of  the  three  lowest  energy  hyperfine  states  of 
6Li  are  stable  against  two-body  loss  processes  for  almost 
the  entire  range  of  field  values  between  15  and  953  G, 
making  this  three-state  mixture  well  suited  to  studies  of 
three-body  physics.  We  observe  a  narrow  loss  feature  at 
127  G  [32]  which  may  be  due  to  a  three-body  resonance 
and  an  additional  state-dependent  loss  feature  at  504  G. 
Finally,  we  measure  the  three-body  loss  rate  coefficient  at 
high  fields  where  all  three  scattering  lengths  are  asymp- 
toting  to  the  triplet  scattering  length  at  =  —  2140ao-  This 
measurement  and  an  observed  trend  toward  lower  loss 
at  higher  fields  has  important  implications  for  realizing 
cold  atom  analogs  of  color  superfluid  and  baryon  phases 
in  QCD. 

We  study  a  6Li  Fermi  gas  with  equal  populations  in  the 
three  lowest  energy  hyperfine  states.  Figure  1(a)  shows 
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FIG.  1:  (color  online)  (a)  Energies  of  the  hyperfine  ground 
states  of  "Li  in  a  magnetic  field.  We  create  a  balanced  mixture 
of  the  three  lowest  energy  states  by  applying  a  RF  pulse  at 
a  Held  of  568  G  (see  inset),  (b)  The  three  pairwise  s-wave 
scattering  lengths  as  a  function  of  magnetic  field  [33]  show 
Feshbach  resonances  at  690,  811,  and  834  G  and  all  asymptote 
to  a,  =  —  2140«o  at  high  field. 


the  energies  of  the  hyperfine  ground  states  of  6Li  as  a 
function  of  magnetic  field.  At  zero  field,  the  three  states 
correspond  to  |1)  =  |A,+A),  |2)  =  |i,-A)  and  |3)  = 
|§,— §)  in  the  | /,  m/)  basis.  For  fields  above  ~  200 G, 
these  states  become  increasingly  electron-spin  polarized 
and  primarily  differ  by  their  nuclear  spin  projection  (i.e. 
in  the  |ms,Tni)  basis  )1)  ~  |  -  ±,1),  |2)  ~  |  -  £,0) 
and  |3)  ^  |  —  5,  — 1)).  The  pairwise  s-wave  scattering 
lengths  between  these  states  (ai2,  <123,  and  <113)  exhibit 
three  broad  overlapping  Feshbach  resonances  as  shown  in 
Fig.  1(b)  [33].  For  very  large  magnetic  fields,  each  of  the 
pairwise  scattering  lengths  asymptote  to  the  triplet  scat¬ 
tering  length  for  6Li,  at  =  — 2140ao.  Therefore,  in  the 
limit  of  large  magnetic  fields,  this  three-state  mixture 
becomes  invariant  under  global  SU(3)  transformations  of 
the  spinor  wavefunction. 

Two-body  spin-flip  processes  are  expected  to  be  espe¬ 
cially  small  for  this  ultracold  6Li  mixture.  Spin-exchange 
collisions  become  energetically  forbidden  even  in  a  small 
magnetic  field.  Dipolar  relaxation  is  small  since  (1)  the 
energy  released  in  a  spin-flip  is  not  enough  to  overcome 
the  centrifugal  barrier  in  the  d- wave  exit  channel  and  (2) 
electron  spin-flip  processes  are  suppressed  at  high  field 
as  the  gas  becomes  increasingly  electron-spin  polarized 
in  the  lowest  energy  electron  spin  state. 

We  can  produce  a  degenerate  Fermi  gas  (DFG)  with 
an  equal  mixture  of  atoms  in  states  |1)  and  |2)  once  ev¬ 
ery  5  seconds  by  evaporatively  cooling  this  mixture  in 


a  red-detuned  optical  dipole  trap  [34].  During  the  first 
second,  ~  108  6Li  atoms  from  a  Zeeman-slowed  atomic 
beam  are  collected  in  a  magneto-optical  trap  (MOT). 
A  crossed  optical  dipole  trap  which  overlaps  the  MOT 
is  then  turned  on  and  optimally  loaded  by  tuning  the 
MOT  laser  beams  ~  6  MHz  below  resonance  and  reduc¬ 
ing  their  intensity  for  7  ms.  The  atoms  are  then  optically 
pumped  into  states  |1)  and  |2)  prior  to  extinguishing  the 
MOT  laser  beams  and  field  gradient.  The  two  beams 
which  form  the  optical  trap  derive  from  a  linearly  polar¬ 
ized  multi-longitudinal-mode  110  W  fiber  laser  operating 
at  1064  nm.  The  beams  nearly  co-propagate  in  the  verti¬ 
cal  direction  ( y ),  have  orthogonal  polarizations  and  cross 
at  an  angle  s;  11°.  The  beams  are  elliptical  with  calcu¬ 
lated  e-2  waist  radii  ~  30pm  and  ~  100pm  at  the  point 
of  intersection.  The  maximum  trap  depth  per  beam  is 
~  1  mK  allowing  ~  5  x  106  atoms  to  be  initially  loaded. 
We  apply  a  bias  field  and  a  noisy  RF  pulse  to  create  a 
50-50  mixture  of  atoms  in  states  |1)  and  |2)  [34]. 

Forced  evaporation  of  the  atoms  occurs  at  a  field 
of  330  G  where  the  two-state  scattering  length  ^ 
— 280ao-  We  lower  the  depth  of  the  optical  trap  Uo  by 
a  factor  of  107  over  3.6  seconds  to  obtain  our  final  tem¬ 
perature.  For  the  work  presented  here  we  study  a  gas 
at  a  temperature  T  >  0.5  7>  (Tp  is  the  Fermi  temper¬ 
ature)  so  that  it  is  appropriate  to  treat  the  cloud  as  a 
thermal  gas.  To  suppress  further  loss  by  evaporation  for 
the  remainder  of  the  experiment,  Uq  is  then  increased  by 
a  factor  of  4  increasing  the  ratio  of  Uo/kpT  by  a  fac¬ 
tor  ~  2.  The  final  oscillation  frequencies  of  the  trap  are 
ux  =  3.84  kHz,  vy  =  106  Hz  and  vz  =  965  Hz  [35]  with  a 
final  trap  depth  per  beam  ~  40  pK.  The  total  number 
of  atoms  ~  3.6  x  105  in  a  balanced  two-state  mixture 
at  T  ^  1.9  pK  for  which  Tp  ~  3.7pK.  The  background 
limited  1/e-lifetime  ~  30  s. 

To  create  an  incoherent  three-state  mixture  we  first 
increase  the  strength  of  the  magnetic  field  in  10  ms  to 
568  G  where  the  mixture  is  stable.  We  then  apply  a  noisy 
RF  magnetic  field  with  two  frequencies  centered  on  the 
|1)  — 12)  and  |2)  — 13)  spin-flip  transitions,  both  broadened 
to  a  width  of  1  MHz  by  frequency  modulation  with  white 
noise.  A  magnetic  field  gradient  of  1  G/cm  in  the 
direction  of  the  bias  field  (z  direction)  is  simultaneously 
applied  to  dephase  the  internal  state  coherence  of  the 
sample.  The  noisy  RF  field  remains  on  for  50  ms  to 
create  an  equal  mixture  of  N  ~  1.2  x  105  atoms  in  each 
state.  At  this  point,  T  =  1.9  pK  with  Tp  =  3.2  pK. 

In  our  first  experiment  to  explore  the  field  dependence 
of  the  loss  mechanisms  in  this  three-state  Fermi  gas,  we 
measured  the  fraction  of  trapped  atoms  remaining  af¬ 
ter  an  evolution  time  of  200  ms  at  fixed  magnetic  fields 
between  15  and  953  G.  Figure  2(a)  shows  the  fraction 
remaining  for  each  of  the  spin  states.  For  each  measure¬ 
ment,  the  magnetic  field  was  ramped  in  10  ms  from  568  G 
to  the  field  of  interest  where  it  was  held  constant  for  ei¬ 
ther  1  or  201  ms.  The  field  was  then  ramped  in  10  ms  to 
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FIG.  2:  (color)  (a)  The  fraction  of  atoms  remaining  in  each  of 
the  three  spin  states  after  200  ms  of  evolution  at  the  magnetic 
field  of  interest,  (b)  Magnetic  field  dependence  of  decay  in 
individual  two-state  mixtures.  (•)  indicates  the  fraction 
remaining  in  state  |2)  (|3>)  of  a  |2)-|3)  (|1)-|3>)  mixture  after 
200  ms.  The  thick  lines  guide  the  eye.  The  representative 
error  bars  indicate  the  standard  deviation  in  the  mean. 


953  G  and  held  for  20  ms  before  a  300  /is  time-of-flight  ex¬ 
pansion  and  absorption  imaging  at  this  field.  The  10  ms 
field  ramp  to  953  G  (where  ai2,  <123, 013  <  0)  ensured  that 
any  atoms  which  had  formed  weakly  bound  molecules  in 
the  vibrational  state  associated  with  the  Feshbach  res¬ 
onances  (but  remained  trapped)  would  be  disassociated 
and  measured.  Further,  each  of  the  spin  states  could  be 
spectroscopically  resolved  at  this  field. 

To  confirm  that  two-body  loss  does  not  contribute  sig¬ 
nificantly  to  the  measurements  above,  we  measured  the 
fraction  of  atoms  remaining  after  an  evolution  time  of  200 
ms  at  fixed  magnetic  field  values  for  each  of  the  three  pos¬ 
sible  two-state  mixtures.  Figure  2(b)  shows  the  fraction 
remaining  in  state  |3)  for  a  |1)  —  |3>  mixture  and  state 
|2)  for  a  |2)  —  |3)  mixture.  The  fraction  of  atoms  remain¬ 
ing  in  the  |1)  —  |2)  mixture  (not  shown)  is  similar  to  the 
|2)  —  |3)  data.  With  the  exception  of  loss  features  be¬ 
tween  600  and  750  G,  each  mixture  was  stable.  The  loss 
features  are  due  to  three-body  recombination  to  bound 
molecular  states  associated  with  the  Feshbach  resonance. 
Since  the  two-state  mixtures  are  stable  for  fields  <  600  G 
and  >  750  G,  we  can  interpret  loss  features  in  a  three- 
state  mixture  at  these  fields  as  being  due  to  three-body 
events  involving  one  atom  from  each  spin  state. 


In  Fig.  2(a),  the  broad  dominant  loss  feature  centered 
at  720  G  occurs  in  the  vicinity  of  three  overlapping  in¬ 
terspecies  Feshbach  resonances.  Significant  loss  due  to 
three-body  recombination  is  expected  near  these  reso¬ 
nances  since  recombination  events  involving  one  atom 
in  each  spin  state  are  not  suppressed  by  the  exclusion 
principle  and  one  expects  a  significant  increase  in  the 
event  rate  when  two  or  more  scattering  lengths  are  res¬ 
onantly  enhanced  [36].  Similarly,  high  stability  at  zero 
field  and  near  the  zero  crossings  of  the  Feshbach  reso¬ 
nances  is  not  surprising  since  at  least  two  of  the  scatter¬ 
ing  lengths  are  small  at  these  fields.  At  high  field  where 
di2,a,23,ai3  — >  at,  the  stability  increases  relative  to  that 
at  720  G. 

Unexpected  resonant  loss  features  are  observed  at  127 
and  504  G,  where  the  two-body  scattering  lengths  are  not 
predicted  to  exhibit  any  resonances.  A  possible  explana¬ 
tion  for  the  feature  at  127  G  is  that  the  binding  energy 
of  a  three-body  bound  state  crosses  through  zero  near 
this  field  value.  At  504  G,  we  observe  differing  degrees  of 
enhanced  loss  for  the  three  states  leading  to  population 
imbalance.  The  feature  at  228  G  in  Fig.  2  is  consistent 
with  our  estimated  location  of  a  |1)  — 13)  p- wave  Feshbach 
resonance.  We  also  observed  a  very  narrow  but  inconsis¬ 
tent  loss  feature  at  259  G  (not  shown). 

A  second  set  of  experiments  measured  three-body  loss 
coefficients  at  fixed  magnetic  fields  where  the  decaying 
populations  remained  balanced.  In  these  experiments, 
the  magnetic  field  is  swept  to  the  field  of  interest  B 
in  10  ms  after  creation  of  the  three-state  mixture  and 
the  number  N(t)  and  temperature  T(t)  of  atoms  in  state 
|3)  are  then  measured  by  time-of-flight  absorption  imag¬ 
ing  at  the  field  B  for  various  delay  times  t.  To  extract 
the  three-body  loss  coefficients,  we  fit  this  data  assuming 
one-  and  three-body  loss  and  heating  due  to  three-body 
recombination. 

When  three-body  recombination  involves  one  atom 
from  each  spin  state,  the  number  of  trapped  atoms  in 
each  of  the  equally  populated  spin  states  N(t)  evolves 
according  to  N  =  —L\N  —  L3(n2)N.  Here,  (n2)  is  the 
average  value  of  the  squared  density  per  spin  state  and 
L3  (Li)  is  the  three-body  (one-body)  atom-loss  rate  co¬ 
efficient.  The  density  distribution  in  the  trap  is  well  de¬ 
scribed  by  a  thermal  distribution  so  that 


dN_ 

dt 


-LxN 


(1) 


where  7  =  L3(mai2/27rfcfl)3/\/27  and  u>  =  (wxu>yu)z)l/3 . 
Following  Ref.  [37],  we  model  the  temperature  increase  in 
the  gas  as  arising  from  “anti-evaporation”  (the  preferen¬ 
tial  loss  of  the  highest  density  atoms)  and  recombination 
heating  7),  which  causes  the  temperature  of  the  atoms  to 
evolve  as 


dT  N2  (T  +  Th) 
dt  ~  7T3  3 
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FIG.  3:  (a)  Three-body  atom-loss  rate  coefficients,  L3,  at  var¬ 
ious  magnetic  field  values.  Evolution  of  (b)  atom  number  and 
(c)  temperature  at  a  particular  field  (180  G).  The  solid  lines 
show  a  best-fit  to  the  experimental  data  which  determines  L3. 


For  a  given  set  of  parameters  (N(0),  T(0),  Li,  L3  and 
T),),  these  two  coupled  differential  equations  can  be  nu¬ 
merically  integrated  to  give  N(t)  and  T(t).  The  best-fit 
parameters  for  a  data  set  at  a  given  field  are  determined 
by  minimizing  x2  with  respect  to  these  fit  parameters. 
For  example,  Fig.  3  shows  the  measured  number  (b)  and 
temperature  (c)  as  a  function  of  time  at  the  field  180  G 
with  the  best-fit  curves.  The  best-fit  L3  values  at  various 
fields  are  displayed  in  Fig.  3(a). 

As  shown  in  Fig.  3(a),  L3  varies  by  over  four  orders 
of  magnitude  and  corroborates  the  features  observed  in 
Fig.  2(a).  The  error  bars  (typically  comparable  to  the 
circles)  indicate  the  uncertainty  in  the  determination  of 
L3  due  to  statistical  fluctuations  and  do  not  include  a 
systematic  fractional  uncertainty  in  L3  of  ±70%  arising 
from  our  fractional  uncertainty  in  atom  number  (±30%) 
and  trap  frequencies  (±5%).  In  the  range  600  G  <  B  < 
750  G,  the  L3  values  (o)  may  be  overestimates  since  the 
measured  decay  includes  loss  events  we  observe  in  two- 
state  mixtures  at  these  field  values.  Near  0  and  568  G, 
three-body  recombination  rates  are  relatively  small  and 
very  stable  three-state  Fermi  gases  can  be  created.  For 
L3  =  10-25  cm6/s,  a  gas  with  a  density  n  =  5  x  1011  cm-3 
per  spin  state  has  a  calculated  lifetime  ~  40  seconds. 
Consistent  with  the  data  shown  in  Fig.  2(a),  a  resonant 
peak  in  L3  is  observed  at  B  ~  127  G.  The  highest  L3 
value  we  report  (10_21cm6/s)  occurs  at  775  G  and  is 
likely  unitarity  limited  (at  these  fields  T  ~  6  pK  due  to 
heating  and  y/ h2/mkBT  =  2200ao).  We  observe  that  L3 
decreases  by  a  factor  of  ~  2.5  as  the  field  is  increased 
from  775  to  953  G,  the  highest  field  we  currently  access. 
This  trend  suggests  that  L3  may  decrease  further  in  the 
high  field  limit  where  ai2,  023,  ai3  — >  -2140ao. 

In  conclusion,  we  can  create  long-lived  three-state 
DFGs  of  6Li  atoms  over  a  range  of  fields.  Additionally, 
the  minimal  two-body  loss  we  observe  makes  this  sys¬ 
tem  promising  for  future  studies  of  Efimov  trimers  and 
other  three-body  bound  states.  Indeed,  the  loss  reso¬ 
nance  at  127  G  may  be  a  signature  of  such  a  state.  The 


recombination  rates  at  high  fields  need  not  preclude  ex¬ 
periments  designed  to  study  color  superfluidity  and  trion 
formation  [18].  For  instance,  by  working  at  low  density 
(n  ~  5  x  1010  cm-3)  and  using  large  period  (d  ~  2  pm) 
optical  lattices,  long  lifetimes  (>  0.1s)  and  strong  inter¬ 
actions  can  simultaneously  be  achieved. 
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please  address  your  inquiries  to  the  Director  of  Contracts.  Should  any  questions  arise  that 
concern  the  scientific  aspects  of  the  research,  please  address  your  inquiries  to  Dr  Anne 
Matsuura,  telephone  703-696-6204. 

Thank  you  for  your  interest  in  the  research  program  of  the  Air  Force  Office  of 
Scientific  Research. 


Sincerely 


Patrick  G.  Carrick 
Director 

Physics  and  Electronics 


The  Basic  Research  Manager  of  the  Air  Force  Research  Laboratory 


